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Truncated Determinant Approach to Light Dynamical Quarks 
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Virtual quark effects in lattice QCD are studied in an approach where the infrared modes are included in a 
precise and gauge-invariant manner. At fixed physical volume the required number of modes does not increase 
in the continuum limit, and the acceptance is tolerable even for very light quark masses. Topological charge 
distributions at varying quark mass, and the static quark potential, showing the onset of string-breaking, are 
presented in both QED2 and QCD4. 
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1. Role of Infrared Dirac Modes in Non- 
perturbative Lattice QCD 

It is well known that the quenched approxima- 
tion seriously mutilates the infrared Dirac struc- 
ture of lattice QCD. Among the many undesirable 
effects of this truncation of the theory one finds 

(1) Nonintegrable quark propagator singulariities 
render the Wilson-Dirac quenched path integral 
undefined [Q, a pathology which surfaces in ex- 
plicit Monte Carlo simulations in the guise of "ex- 
ceptional configurations" leading to uncontrolled 
statistical fluctuations. 

(2) Quenched chiral logarithms complicate the ex- 
traction of reliable chiral extrapolations. 

(3) Nontrivial topology (i.e. nonzero topological 
charge) is not suppressed in the light quark limit. 

(4) Characteristic physical phenomena induced 
by virtual quark pairs- string-breaking, rho de- 
cay, etc- are absent. 

Inclusion of the infrared quark modes, character- 
ized as the low eigenvalues of the hermitian opera- 
tor H = 75 (Ijt — m), allows a precise gauge invari- 
ant truncation which interpolates in a natural 
way between the quenched and fully unquenched 
theory, eliminating all the aforesaid problems. 
The spectrum of H can be thought of as char- 
acterizing the gauge-invariant off-shellness of the 
quark field, as the eigenvalues are gauge-invariant 
and reduce in the free theory to ±\/p 2 + m 2 for 
a quark mode of Euclidean momentum p. We 



shall discuss simulations in which all modes up to 
a cutoff somewhat greater than (typically, about 
twice) Aqcd are included in the quark determi- 
nant part of the full Boltzmann weight. It turns 
out that accurate inclusion of the remaining ul- 
traviolet modes ||] requires a cutoff at a fixed 
number of eigenvalues rather than in energy. 

2. Implementing a Gauge Invariant Trun- 
cation of the Quark Determinant 

The spectrum of the hermitian Dirac operator 
H consists of eigenvalues Xi(A) which are gauge- 
invariant functions of the gauge field A. For 
QCD4 the dimension of H is N = 12V, with 
V the lattice volume, and it is convenient to or- 
der the spectrum with the index i running from 
—Af/2 to Af/2 with Aj < A i+ i. With the conven- 
tional lattice normalization the spectrum extends 
from -(1 + 8k) to +(1 + 8k) (roughly, from -2 
to +2) with the physical branch extending up to 
|A| ~ 0.5. We shall describe the results of simula- 
tions in which the full quark determinant contri- 
bution is replaced by T>n x (A) = Ylt=^ N \(A). 
For QCD4 we choose N\ ~ 50, corresponding to 
An x ~ 370 Mev, in other words, somewhat larger 
than the low energy scale Aqcd- Low eigenvalues 
are obtained by the following Lanczos scheme: 
(1) Starting with a randomly chosen initial vec- 
tor v\, an orthonormal sequence v\,V2,v^, ..vn l 
is generated by the standard recursion: 
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Figure 1. Pseudo-scalar correlators in QED2 



v k+1 = j- k Hv k - ^« fc -i - f k v k 
with the constants a k , f] k determined from over- 
laps of generated vectors. In the basis of the Vi, H 
is tridiagonal. The corresponding real symmetric 
tridiagonal matrix Tn l has the a k on the diago- 
nal and the (5 k on the sub (and super) diagonal. 

(2) A Cullum-Willoughby sieve j| is applied to 
remove spurious eigenvalues. 

(3) The remaining "good" eigenvalues converge 
most rapidly in the needed infrared portion of the 
spectrum. The stability and accuracy of the con- 
verged eigenvalues has been checked extensively 
by gauge transforming the gauge field. 

(4) The diagonalization of the Tm l matrix (typ- 
ically, of order 10,000 in the QCD case) can be 
completely parallelized using the Sturm sequence 
property Q of tridiagonal matrices which allows 
nonovcrlapping parts of the spectrum to be inde- 
pendently extracted by a bisection procedure. 

3. Truncated Determinant Simulations in 
QED2 

We have tested the truncated determinant ap- 
proach in two-dimensional abelian gauge theory 
(QED2) on 10x10 and 16x16 lattices at (3 =4.5 
and for bare quark masses mo =0.06 and 0.10. 
Superrenormalizability implies that higher mo- 
mentum modes are essentially inert (the prod- 
uct of the upper 90% of the spectrum is prac- 
tically constant) and one finds complete agree- 
ment between exact dynamical simulations and 
those using only the 2A^ =20 lowest eigenmodes 
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Figure 2. Qt op histogram in QED2 



Quark-antiquark potentials (QED2) 
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Figure 3. Quark-antiquark potentials in QED2 



of the (hermitian) Wilson-Dirac operator. The 
simulation procedure used was as follows. A pure 
gauge update consisting of 5 metropolis sweeps 
was followed by an accept-reject step using the 
recomputed value of the truncated determinant 
T>m x (A). Strictly speaking || detailed balance 
only holds for the pure gauge step if links or sets 
of noninterfcring links are updated in random or- 
der. For our QED2 runs we have found no statis- 
tically significant differences between the results 
obtained with ordered or random link updates. 
(A fully parallel code incorporating detailed bal- 
ance will be used in all future truncated determi- 
nant simulations). The acceptance ratio for this 
procedure is tolerable (typically from 20-40%) be- 
cause the typical fluctuation in logT>N x {A) after 
a metropolis sweep is of order unity. (This is also 
the case in QCD4 with N\ chosen large enough 
to include all the low energy chiral physics). 
Some results obtained from truncated determi- 
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Figure 4. Q top histogram in QCD4 for k=0.1597 
(diamonds- chiral result) 
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Figure 5. Static potential for quenched and dy- 
namical QCD4. Error bars are largest for lightest 
mass and only shown for this case 



nant simulations in QED2 arc shown in Figs. 1- 

3. In Fig(l) the "pion" (i.e. pseudoscalar) cor- 
relators are shown for a very light quark: inclu- 
sion of the lowest 20 eigenvalues leads to correla- 
tors essentially identical to those in the fully un- 
quenched theory. In Fig(2) the topological charge 
(Qto P ) distribution is shown for the same param- 
eters in quenched, unquenched and truncated de- 
terminant versions of the theory. The suppression 
of nonzero topological charge is basically com- 
plete when only 5% of the spectrum is included in 
the determinant. Fig. 3 shows the onset of string 
breaking in the quark-antiquark potential mea- 
sured from Wilson loops on a 16x16 lattice. 

4. Truncated Determinant Simulations in 
QCD4 

We have used the truncated determinant ap- 
proach in QCD4 on a 12 3 x24 lattice at ,3=5.9 
for three kappa values, k=0.1570, 0.1587 and 
0.1597. Pion correlators were measured in order 
to perform a chiral extrapolation and determine 
the critical k c . 100 modes (up to ~370 MeV) 
were included in the quark determinant, with ei- 
ther one or two complete Metropolis pure gauge 
sweeps between determinant accept-reject steps 
(with acceptance varying betwen 20% and 50%). 
The topological charge can be extracted from an 
anomalous chiral identity, and is related to the 
trace of the inverse of H = 75 (Lj) — m) . This spec- 
tral sum converges rapidly: summing the lowest 



100 eigenvalues measures the topological charge 
(as per this definition) to a few percent. The 
histogram of measured topological charge in the 
lightest case, k=0.1597 (corresponding to a pion 
mass of about 280 MeV) is shown in Fig. 4, to- 
gether with the analytic prediction for this distri- 
bution following from a chiral analysis [Q. 

In Fig. 5 the static quark-antiquark potential 
(obtained from the correlation of Wilson lines in 
Coulomb gauge) shows a clear screening effect 
although true asymptotic flattening presumably 
occurs at distances where statistical fluctuations 
and finite volume effects dominate. 

REFERENCES 



W. Bardee n, A. Duncan, 



E. Eichten and H. 

Thacker, |hep-lat/9806002 . 
Talk of A. Duncan, Lattice 97, Proc. Suppl. 
Nucl. Phys. B63, 811 (1 998); A. Duncan E. 



Eichten and H. Thacker, frep-lat 9806"020 



Talk of E. Eichten, this conference. 
J. Cullum and R.A. Willoughby, J. Comp. 
Phys. 44, 329 (1981). 

G. H. Golub and C.F. Loan, Matrix Compu- 
tations, 2nd edition (Johns Hopkins, 1990). 
We thank D. Toussaint for explaining the con- 
nection of the detailed balance condition to 
the link update procedure. 

H. Leutwyler and A. Smilga, Phys. Rev. D46, 
5607 (1992). 



